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We propose a correlated spin-singlet-pairs wave function to describe the spin-gap phase of the 
one-dimensional t — J model at low density. Adding a Jastrow factor with a variational parameter, 
u, first introduced by Hellberg and Mele M, is shown to correctly describe the long-range behavior 
expected for the Luther-Emery phase. Using the variational Monte Carlo method we establish a 
relation between v and the Luttinger exponent Kp, f 



j_ 



74.20.-Z, 75.10.Jm, 75.50.Ee 



0^ 
O 



X3 

o 



X 



The high temperature superconductivity(HTSC) has 
continued to be one of the important current issues in 
the field of condensed matter physics. It is generally be- 
lieved that the two-dimensional (2D) t—J model is a good 
working model for the theory of HTSC, which involves 
the basic interactions of the copper-oxygen planes. How- 
ever, recently the one-dimensional (ID) version of the 
model has received much attentions. Inspired by the un- 
usual normal-state properties of high temperature super- 
conductors (HTS), Anderson ||^ proposed that HTS may 
be described as the Tomonaga-Luttinger liquid (TLL) ||] 
instead of the conventional Landau's Fermi liquid. Since 
the TLL is well studied in the ID model, many of its 
properties give us useful references in the study of the 
2D model. For example, the shift of the characteristic 
momentum from 2kp to 2kp^ of spinless fermions (SF) 
in the density-density correlation function has been used 
to argue for spin and charge separation |^,^. Recently, 
there are more experimental evidences for a spin-gap 
phase Q in HTS. So far there is no good account of this 
phase from the 2D model. It turns out that there is also 
a phase with a spin gap in the ID model. Recent stud- 
ies have identified this phase as the Luther-Emery 
(LE) phase A more careful study of this phase in ID 
would help us to gain insights in dealing with the 2D. In 
this paper we will present a wave function that catches 
the essence of this LE phase. 

The ID t ~ J model is defined as 



Htj = -t^(c+Ci+iCT + h.c.) 



(1) 



with the constraint of no double occupancy. By diag- 
onalizing the 16-site chain, Ogata et al. |]l^ has found 
two phases in the phase diagram of the ID t — J model. 
For very large J/t holes and spins phase separate. With 
decreasing J the system is described as a TLL which ex- 
hibits power-law correlations with exponents charater- 
ized by a single parameter, Kp. For TLL there is no gap 



in both spin and charge excitations. In addition to the 
momentum distribution, the important correlation func- 
tions such as spin-spin, density-density and pair-pair, all 
have power-law decay: 

{S,{r)SM) - Aor-^ + Aicos{2kFr)r~^' 

+ A2Cos{4:kFr)r~^\ (2) 
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and some of the ex- 



where A(i) = Ci^ Ci- 
ponents are: Xs — Kp + 1 , Xj\i = Kp + 1 and Ap = + 1 . 
Using the ground-state projection method, two groups 
have recently found a third phase, the Luther- 
Emery(LE) |9j phase with a non-zero spin gap in the re- 
gion of low electron density and high interaction strength. 
Due to the spin gap, the spin-spin correlation function 
decays exponentially with distance. It has similar power- 
law correlation functions for density-density and pair- 
pair as the TLL but with different exponents: A^v = Kp 
and Ap = -A-. 

Variational approaches have been very successful in the 
study of the phase diagram. In particular, Hellberg and 
Mele (HM) H have shown that the TLL is very well 
represented by a wave function with long-range Jastrow 
correlations. The HM wave function is defined as 



\HM) = H, 
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where denotes the hole position, L is the total number 
of sites and $p is the projected ideal Fermi gas wave 
function. The holes repel each other when v is positive 
and attract otherwise. HM also showed that with tuning 
u for different J/t and densities this wave function can 
reproduce the phase diagram of the IT) t — J model. Not 
only the energy is quite accurate it also produces the 
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correct power-law correlations that are the signatures of 
the TLL. Specifically they have found 
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However \HM) fails to predict the spin-gap phase 
which is a LE phase instead of a TLL. 

In the study of the spin-gap phase, using the exact 
pairing wave function of two electrons in an infinite chain 
as a basis, Chen and Lee (CL) proposed a singlet-pair 
(SP) wave function. 



(7) 



where h = 2t/J. Ne is the total number of electrons, 
Pd is the projection operator that forbids two parti- 
cles occupying the same site and the operator = 
Y.i'^ii'^i+ni ~ ^ii^i+nf Without any tuning param- 
eters CL showed that \SP) has lower energy than \HM) 
and more significantly it has the correct short distance 
spin-spin correlation which characterizes the LE phase. 
However with increasing particle density there will be 
correlations among pairs and holes. And more seriously, 
\SP) is of a particular form of the projected BCS wave 
function (or RVB) hence produces a long-range pairing 
order, which is inconsistent with the LE phase with a 
power-law correlation in pairing. Taking h as a tuning 
parameter will not help to suppress the long-range order. 

In view of the power-law correlation produced by the 
Jastrow factor in \HM), a natural way to modify the 
wave function \SP) in order to correctly represent the 
pairing correlation in the LE phase is to add the Jastrow 
factor to \SP). This new trial function denoted as \CSP) 
is 



\CSP) = n,: 
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\SP). 



(8) 



Here we present our results of \CSP) by using the vari- 
ational Monte Carlo method. The closed-shell boundary 
condition is used for all the data presented in this pa- 
per. The correlation functions cither in real space or 
in Fourier space are shown in Fig.l for various v. The 
pairing correlation function is defined in Eq.(4). The 
variational parameter is ft, = 2/3 for J/t ~ 3 at the elec- 
tron density = 1/6. Open circles are the results of 
\HM) with V ~ —0.5 which is the optimal wave func- 
tion within \HM) in the spin-gap region of ID t — J 
model. Comparing the variational energies among sev- 
eral wave functions for J/t = 3 we find that the lowest 
energy is obtained by \CSP) with u — Q.l (open trian- 
gles). Fig. 1(a) shows that in this region of J/t, \HM) 
underestimates the pairing correlation as compared to 
the \CSP). This has been pointed out in previous stud- 
ies using the ground-state projection method |0,D. On 



the other hand, the long-distance bahavior of the pairing 
correlation of \SP) {v = 0, filled circle) seems to lead to 
a long-range pair ordering as mentioned above. With a 
non-zero u the long-range behavior has changed and we 
will show below that it has the power-law dependence as 
expected from the LE phase. 
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FIG. 1. Variational Monte Carlo results of (a)pairing cor- 
relations (shown in real space), (b)spin and (c)density struc- 
ture factors of \HM,v — —0.5) and \CSP) for several differ- 
ent u. The system has 10 electrons in 60 sites. Symbols are 
defined in Fig. 1(b). 

Spin {S{q)) and density {N{q)) structure factors are 
plotted in Figs. 1(b) and 1(c) respectively, with the same 
parameters as in Fig. 1(a). It is easy to show that S{q) 
should be quadratic at small q if there is a gap in spin 
excitations [0. This indeed has been found in Fig. 1(b) 
for \SP) and \CSP). More interestingly, this behavior is 
quite robust, it hardly changes with the variation of the 
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exponent v in contrast to the dramastic changes occurred 
in pairing and density correlations. The reason, we be- 
Ueve, is that the wave functions \SP) and \CSP) only 
have short-range pairs. This preserves the spin gap even 
in the presence of strong density fluctuations induced by 
the Jastrow correlation factor. 
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FIG. 2. Finite-size scaling of the quantities (a)6p(I//2) 
and (b)6jv(L/2) as defined in Eqs.(9) and (10) for \CSP) with 
a fixed h = 0.6. The electronic density is fixed at rie = | for 
lattice sizes ranging from L =36 to 132. The solid lines are 
the least-square fits of the data points. The slops of the lines 
in (a) and (b) give 1 — Ap and 1 — Ajv, respectively, (c) is 
similar to (a) but for a fixed v = 0.2 and various h. 

As shown in Fig. 1(c) |7JM, = —0.5) has the maxi- 
mum of iV((7) at 2-k/L. This indicates the system is close 
to the phase separation. However for the LE phase one 
would expect a sharp peak at 2kF and even a divergent 
one if Kp < 1. For \SP) the maximum is indeed at 2kp 
but it is quite broad. As for \CSP) we find that the peak 



becomes sharper and sharper as u increases. The size of 
the peak also grows with u. By using finite-size scaling 
we find that the peak diverges with the lattice size for 
v > 0.5 (see below). 

In addition to establish the power-law behavior for 
density-density and pair-pair correlations we must also 
find a relation between their exponents. In order to 
extract the exponents of the power-law correlations we 
consider the finite-size scaling behavior of the following 
quantities |12[, 



bp{L/2)^P{q = 0)-P{q = 2TT/L) (9) 



and 



bN{L/2) ^ 2N{q ^ 2kF) - N{q = 2kF + 2it / L) 

- N{q^2kF - 2tt/L), (10) 

where P{q) is the Fourier transform of P{r). In Fig. 2(a) 
we plot log \bp{q = 0, L/2)] versus log(L/2) to obtain the 
exponent 1 — Ap as the slope of the linear fit of the data. 
A similar analysis oibN^q = 2kF, L/2) gives the exponent 
1 — Aat. We have used lattice sizes ranging from L=36 to 
132. The very well linear fit of these plots supports our 
conclusion of power-law behavior in pairing and density 
correlations. The exponents obtained in these plots are 
tabulated in Table I. We find that the relation Ap = is 
satisfied with the data in Table I, i.e., the scaling relation 
expected for the LE phase is recovered Q . Additionally 
we find that within the error bars the variational param- 
eter 1/ is related to the critical exponent Kp of the LE 
phase in a simple relation p3| 



2i^ 



(11) 



Correspondingly we have Xp = and A^v = 2u 



TABLE I. Critical exponents Ap(AAr) of pairing (den- 
sity) correlations for \CSP) with various varational pa- 
rameters ly. The numbers in the parentheses show the 
error in the last-digit. The last row is the predicted val- 
ues oiKp by Eq.(ll). 
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0.21(2) 0.43(3) .62(2) 0.77(8) 1.03(4) 1.16(8) 



2.3(4) 1.63(8) 1.32(8) 1.06(8) 0.84(5) 



2.5 1.667 1.25 



1.0 0.833 



The finite-size scaling of pairing correlations for a given 
v{=0.2) and various h is shown in Fig. 2(c). We find that 
all the data fit to lines of the same slope hence of the 
same exponent. Therefore we conclude that the long- 
range power-law correlations are controlled by the Jas- 
trow factors, irrespective of the parameter h which con- 
trolls the short-range properties. 
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FIG. 3. (a) Variational energies, (b)Density and (c) spin 
structure factor, as well as (d) pairing correlation of the t — J 
model for \CSP), \SP) and \HM) at the density = 1/6. 
The variational parameters are optimized for J/t = 2.8. We 
take h = 2t/J for both \CSP) and \SP). And u = 0.08, -0.5 
for \CSP} and \HM), respectively. The lines and symbols are 
labled in (c). This system with 4 electrons in a 24-site chain 
is also exactly diagonalized by the Lanczos method and these 
exact results are shown as the circles. 

Having established \CSP) as a good wave function to 
describe the LE phase, now we like to investigate if it 
is also a proper wave function to faithfully represent the 
ground state. The varational energies of A/), \SP) and 
\CSP) at the density rie = 1/6 are shown in Fig. 3(a). 
We observe that \CSP) and \SP) are very close in en- 
ergy for large J/t. By using the power method we find 
that the variational energy of \CSP) is about 0.3% above 
the ground state energy. Although \CSP) and \SP) have 
similar energies, they have significant differences in their 
long-range correlations as indicated above. In Fig. 3(b) - 
(d) we show the density and spin structure factors and 
pairing correlations for \HM), \SP) and \CSP) opti- 
mized for t — J model with J/t = 2.8. This system 
contains 4 electrons in a 24-site chain. For this system 
the exact ground state can be obtained by the Lanczos 
diagonalization method. Results for the ground state are 
shown as the open circles in Fig. 3. We find that for 
spin correlation \CSP) and \SP) have similar quadratic 
behavior at small q and agree with the exact result very 
well. However for density and pairing correlations \CSP) 
describes much better than \SP). The excellent consis- 
tency between the correlations of the ground state and 
\CSP) supports the expectation of the LE phase in the 
low density region of the ID t — J model. 



In conclusion, we have presented a wave function to de- 
scribe the LE phase in ID. This wave function \CSP) has 
correlated spin-singlet pairs. This is established by using 
the variational Monte Carlo method. The wave function 
shows exponential dependence for the spin-spin correla- 
tion and power-law behavior in density-desity and pair- 
pair correlations. By finite-size scaling we established the 
relation of the varational parameter v to the exponent 
Kp, Kp — Comparing with the exact ground state of 
a small lattice we showed that \CSP) describes very well 
the ground state properties of the spin-gap phase of the 
ID t - J model. 
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